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Abstract 

O - 

£\j , In this article, we study the mass spectrum of the baryon-antibaryon bound states 

pp, EE, HE, AA, pN(14A0), ES(1660), EE' and AA(1600) with the Bethe-Salpeter 
equation. The numerical results indicate that the pp, EE, HE, piV(1440), EE(1660), 
EE' bound states maybe exist, and the new resonances X(1835) and A(2370) can be 
tentatively identified as the pp and p7V(1440) (or N(1400)p) bound states respectively 
^vq ' with some gluon constituents, and the new resonance X(2120) may be a pseudoscalar 

glueball. On the other hand, the Regge trajectory favors identifying the A(1835), 
A(2120) and A(2370) as the excited 7/ (958) mesons with the radial quantum numbers 
n = 3, 4 and 5, respectively. 
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^ 1 Introduction 

in , 

^SJ , In 2003, the BES collaboration observed a significant narrow near-threshold enhancement 

in the proton-antiproton (pp) invariant mass spectrum in the radiative decay J/tp — > jpp 
PP. The enhancement can be fitted with either an S-wave or a P-wave Breit-Wigner 
resonance function. In the case of the S-wave fitted form, the mass and the width are 
M = (I859±f ±| 5 ) MeV and T < 30MeV, respectively. In 2005, the BES collaboration 
observed a resonance state X(1835) in the r]'ir + ir~ invariant mass spectrum in the process 
J/ip -)• jTT + ir-r]' with the Breit-Wigner mass M = (1833.7 ±6.2 ±2.7) MeV and the width 
T = (67.7 ± 20.3 ± 7.7) MeV, respectively [2]. Recently, the Z(1835) was confirmed by 
the BES collaboration in the radiative decay J ftp — > "fir + Tr~ri' with a statistical signifi- 
cance larger than 20 a, the fitted mass and width are M = (1836.5 ± 3.0±|;f) MeV and 
T = (190 ± 9tgg) MeV, respectively 0g|. The mass is consistent with the BESII result 
[2], while the width is significantly larger. Furthermore, the BES collaboration observed 
two new resonances V(2120) and V(2370) in the ir + TT~r] / invariant mass spectrum with 
statistical significances larger than 7.2 a and 6.4 a, respectively. The measured masses 
and widths are M x{2m = (2122.4 ± 6.7^;^) MeV, M x(2370) = (23 76.3 ± 8.7±|;|) MeV, 
Tx(2i20) = (83 ± 16±n)MeV and r x(2 37o) = (83 ± 17+^ 4 ) MeV, respectively [3j S] - 

Many theoretical works were stimulated to interpret the nature and the structure of 
the new resonance X(1835), such as the pp bound state [5], the pseudoscalar glueball [6j[7], 
and the radial excitation of the rj' [HI O [10] , the threshold cusp [H], etc. 

In Ref.[12|. Liu, Ding and Yan study the decay widths of the second and the third 
radial excitations of the pseudoscalar mesons rj and 7/ with the 3 Po model, and observe 
that the interpretation of the 7/(1760) and X(1835) as the second radial excitations of the rj 
and rj' crucially depends on the measured mass and width of the 77(1760), and suggest that 



- 1—1 

X 



1 E-mail, wangzgyiti@yahoo.com.cn. 



1 



there may be sizable pp content in the X(1835), and the AT(2120) and X(2370) cannot 
be understood as the third radial excitations of the r\ and rj respectively, the X(2370) 
probably is a mixture of the r/'(4 1 S'o) and glueball. 

In Ref.[13], we take the AT(1835) as a baryonium state with the quantum numbers 
J PC = h , and calculate the mass spectrum of the baryon-antibaryon bound states pp, 
EE, EE, and AA in the framework of the Bethe-Salpeter equation with a phenomenological 
potential. The numerical results indicate that the pp, EE and EE bound states maybe 
exist, and the X(1835) can be tentatively identified as the pp bound state. In this article, 
we extend our previous work to study whether or not there exist the baryon-antibaryon 
bound states pA>(1440), ££(1660), EH' and AA(1600) (or JV(1440)p, £(1660)£, S'S and 
A(1600)A), here the / denotes the first radial excited state. In the scenario of the coupled 
channel effects or the hadronic dressing mechanism [14] . the pseudoscalar mesons X(1835), 
X(2120) and AT(2370) can be taken as having small pseudoscalar qq kernels of the typical qq 
meson size. The strong couplings to the virtual intermediate hadronic states (for example, 
the pp, piV(1440), EE, EE, etc) can change the bare masses which originate from the 
quark-gluon interactions significantly, and enrich the pure qq states with some baryon- 
antibaryon components. In the present case, the input parameters of the Bethe-Salpeter 
equation should be readjusted to obtain reasonable predicted masses if the coupled-channel 
effects are considered. However, it is beyond the scope of the present work. Such a scenario 
needs detailed studies. 

The article is arranged as follows: we solve the Bethe-Salpeter equation for the baryon- 
antibaryon bound states in Sec. 2; in Sec. 3, we present the numerical results and discus- 
sions; and Sec. 4 is reserved for our conclusions. 



2 Bethe-Salpeter equation 

The Bethe-Salpeter equation is a conventional approach in dealing with the two-body 
relativistic bound state problems, and has given many successful descriptions of the hadron 
properties [15j[16]- We write down the ladder Bethe-Salpeter equation for the pseudoscalar 
bound states in the Euclidean spacetime^l, which can be derived from the Euclidean path- 
integral formulation of the theory, 

/d 4 k 
j^^x(k,Ph 5 G(q-k), (1) 

Sy\(q±Z 1/2 P) = i( 1 -q±t l/2l -P)+M l/2 , 

Mi /2 

^ 1/2 ~ Mi + M 2 ' 

where the P„ is the four-momentum of the center of mass of the baryon-antibaryon bound 
state, the q^ is the relative four-momentum between the baryon and antibaryon, 75 is 
the bare baryon-meson vertex, the x(q> P) is the Bethe-Salpeter amplitude of the baryon- 
antibaryon bound state, and the G{q — k) is the interaction kernel. 

In this article, we use the metric 5 M „ = (1,1,1,1), {7 M 7„ + Jvfn} = 2(5 M „, the momentums fc M = 
{k 4 ,t), g M = (q 4 , V) and P M = (iE,~P) with P 2 = -M%. 
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In the flavor SU (3) symmetry limit, the interactions among the ground state octet 
baryons and the pseudoscalar mesons can be described by the lagrangian C, 



where 



C = V2 (DTr (B {P, B} + ) + FTr (B [P, B}_)) 



B 



P = 



\ 



E+ 



p \ 

n 



V2' 



K~ 



K° 



■TP 



(2) 



(3) 



and the D and F are two parameters for the coupling constants. From the lagrangian, we 
can obtain 



9n°pp — -9n°nn — D + F , # 7r S + S + — —g n O%-x 2F , 

r> r, D - 3F 

9tt°s-e- = -9n s°s = D - F , g VP p = g vnn = 



V3 



3r?E+E+ 



2D 

9riZ Q Z° — ~9vAA — , 

D + 3F 



and write down the kernel G(q — k) explicitly, 



G , q _ k) = dHq-k)^ | g 2 (q-k)C v ^ 
(q — k) 2 + m 2 (q — k) 2 + m 2 ' 



(4) 



(5) 



where the_coefficients C n = (1 + a) 2 , 4a 2 , (1 - a) 2 , and C v = a) " , |, (1+ f* )2 , | for 
the pp, EE, SS, AA bound states respectively; g 2 (k) = D 2 and a = -p. 
With a simple replacement 



B ^ B' , D ^ D' , F ^ F' , a' , g ^ g' 



where 
B' 



( ^E(1660) + TgA(1600) 



E(1660) 



iV(1440) \ 



V 



E(1660) 
£'(?) 



-^E(1660) + ^gA(1600) A^(1440) 

-^A(ieoo) y 



(6) 



,(7) 



we can obtain the corresponding couplings among the first radial excited octet baryons and 
the pseudoscalar mesons. The corresponding coefficients are C n = (1 + a){\ + a'), Acta', 
(l-a)(l-a'), and C v = £=3»Xi=3gQ |, (i+3a)(i+3a') > 4 for the bound stateg pi y(1440), 
EE(1660), HH', AA(1600) (or AT(1440)p, E(1660)E, H'S, A(1600)A), respectively In this 
article, we can take the approximation a = a' and g = g' (i.e. we assume that the ground 
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state and first radial excited octet baryons have the same quantum numbers except for the 
masses), and use the parameters determined in our previous work [T3]. Furthermore, we 
also solve the Bethe-Salpeter equations with the parameters g 1 ^ g, as the flavor SU(3) 
symmetry does not warrant g' ~ g. 

In this article, we have neglected the baryon-antibaryon annihilation effects. There are 
two typical diagrams for the annihilation contributions, the exchanges of the intermediate 
gluons and the intermediate mesons. The contributions of the intermediate gluons oc a™ 
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with n > 2, where a s = f 2 -, the g s is the quark-gluon coupling constant, the values of the a™ 
are usually very small and can be neglected for the heavy-quark systems [17] ■ While in the 
light-quark systems, the gluonic annihilation effects are expected to play a more important 
role. However, even in the light mesons, the annihilation Hamiltonian is usually small, 
and there is considerable phenomenological evidence to reinforce one's expectation that it 
becomes weaker as a meson system becomes more excited [T7j. The effects are somewhat 
alleviated by this factor, thus in practice the gluonic annihilation effects are often ignored. 

The annihilations to the intermediate mesons can be estimated as oc -tjt- ; — r — with 

n > 1. If we take the typical momentum q = y/2M r \E\ ~ 200 MeV, where the M r is the 
reduced mass in the pp system, and the E is the bound energy of the X(1835) as the pp 

bound state, My 3> q, and ( tttt—w^ — ) *C -? 9 — , such annihilation effects can be 
neglected. The annihilation effects in other channels can also be neglected with analogous 
arguments. 

In Ref.|13]. we choose the value a = 0.6 from the analysis of the hyperon semi-leptonic 
decays |18| . and take the coupling constant g 2 (k) as a modified Gaussian distribution 

g 2 (k) = A (^j^j exp m * ne Euclidean spacetime, where the strength A and the 

distribution width fi are two free parameters. The ultraviolet behavior of the modified 
Gaussian distribution warrants the integral in the Bethe-Salpeter equation is convergent. 

The Euclidean Bethe-Salpeter amplitude of the pseudoscalar baryon-antibaryon bound 
state can be decomposed as [16] 



X (q,P) = 1 s{F{q,P)+iPF 1 {q,P)+i4F 2 {q,P) + \4,P]Fz{q,P)} . (8) 

As in our previous work, we take the approximation 

X (q,P) = l5 {F(q,P)+iPF 1 (q,P)} , (9) 

for simplicity. We can denote the baryon fields as ^(x), and perform the Fierz re-ordering 
to study the contributions from different spinor structures^, 

+ \{l^5)a^(x)j fl 7 5 ^(0) + ^(i7 5 ) a i3^(x)i^(0) . (10) 

The pseudoscalar current ^(x)ij^(O), the axialvector current ^(x)'y^ 75^(0) and the 
tensor current ^>(x)a^ u ^f(0) have nonvanishing couplings with the pseudoscalar meson 



3 Here we choose the X(1835) as an example to illustrate the estimation, and use the metric in the 
Minkowski spacetime. 
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X(1835), 



i 75 (0|*(x)i7 5 *(0)|X(P)) oc 75^0 + 759 • -P^! + ■ • ■ , 
7^75(0|^(x) 7At 75^(0)|X(P)) a 75 P5o + 75 + • • • , (H) 
a^(0\y(x)a^(0)\X(P)) oc 75 ^,P]C + ---, (12) 

the coefficients Ao, A\, Bq, B\, Co are functions of the and P„. We can take the 
estimations P M ~ M x ss 1.8 GeV, g M ~ y^M^I « 0.2 GeV, \^,P] ~ y/2M r |£|M x « 
0.36 GeV 2 , translate the P M and into dimensionless quantities, P M — > P^A, — > q^A 
with A = 1 GeV, and absorb the A into the coefficients Aq, A\, Bq, B\, Cq, then P^ ~ 1.8, 



q„ ~ 0.2, 



0.36. Compared with the term 75 P, the term 75 ^ is greatly suppressed 
and can be neglected. Furthermore, we expect that couplings of the tensor currents to the 
pseudoscalar mesons are weaker than that of the pseudoscalar and axialvector currents, 
and neglect the term 75 y, P] , which is also suppressed as the 75 

The Bethe-Salpeter amplitudes F(q,P) and Fi(q,P) can be expanded in terms of 

1 

Tchebychev polynomials T n 2 (cos 0) [19] , 

00 1 
F{q,P) = ^ n F n (g 2 ,PVP n T n 3 (cos0) , 

^ I 

F 1 (q,P) = i n F?(<l 2 , P 2 )q n P n n (cos 0) , (13) 



n=0 



where 6 is the included angle between q^ and P„. If we translate the momenta q^ 
and P M into the dimensionless quantities q^ and P^ respectively, and absorb the A into 

the F n (q 2 ,P 2 ) and Ff(g 2 ,P 2 ), then qPT?(cos0) ~ 0.36 cos ~ 0, c?P 2 T*(cos6) ~ 
0.13 cos 20 ~ —0.03, here we have taken the average cos n0 ~ — cos n0sin 2 0d0. It is im- 
possible to solve an infinite series of coupled equations of the F n (q 2 ,P 2 ) and F™(q 2 ,P 2 ), 
we have to make truncation in one or the other ways. In this article, we neglect the 
small terms with n > 1. Numerical calculations indicate that taking only the terms with 
n = can give satisfactory results. If we take into account the small terms with n > 1, 
the predictions may be improved mildly. In the following, we will smear the index for 
simplicity. 

Multiplying both sides of the Bethe-Salpeter equation by 75 \rf,P] and carrying out the 
trace in the Dirac spinor space, we can obtain an simple relation F = (Mi + M^jFi^ the 
amplitudes F(q 2 ,P 2 ) and F\{q 2 ,P 2 ) are not independent. The Bethe-Salpeter amplitude 
can be written as 

x(q,P) = T 5 (l + 1 ^-)^ 2 ,P 2 ) ) (14) 
and the Bethe-Salpeter equation can be projected into the following form, 

{" 2+M ^ + wrh p2 ) F ^ p2 '> - <i5) 
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We can introduce a parameter X(P 2 ) and solve above equation as an eigenvalue prob- 
lem. If there really exists a bound state in the pseudoscalar channel, the mass of the 
bound state X can be determined by the condition X(P 2 = —M x ) = 1, 

If we take q 2 = and assume that there exists a physical solution, then 

{ M ^-TK^h Mi ) F( °-- Mi) = S§f^\-Ml)G ( 0-k). (17) 

In numerical calculations, we observe that the Bethe-Salpeter amplitude F(k 2 , —M x ) has 
the same sign in the region k 2 > 0, 



MiM 2 n/r2 f d 4 k F(k 2 ,-M 



(Mi + M 2 ) 2 J (2vr) 4 F(0, — M 



f2 > 



M 1 M 2 - — ' = / . o ^ ' AZ G (0-fe)>0, (18) 
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and obtain an simple relation (or constraint), 

M 2 X < (Mi + M 2 ) 2 , (19) 

which survives for q 2 > (although the relation is not explicit for q 2 > 0), i.e. the bound 
energy Ex originates from the interacting kernel G(k) and should be negative, Ex = Mx — 
Mi — M 2 < 0. On the other hand, if the Bethe-Salpeter amplitude F(k 2 , —M x ) changes 

sign in the region k 2 > 0, which does not warrant the positive value J J5^~pn^, m*) G{§ — 

k) > 0, and the relation M 2 X < (Mi + M 2 ) 2 fails to survive. 



3 Numerical results and discussions 

The input parameters are taken as m n = 135 MeV, m v = 548 MeV, M p = 938.3 MeV, 
M s+ = 1189.4 MeV, M s _ = 1321.7 MeV, M A = 1115.7 MeV, M N{1U0) = (1420 - 
1470) MeV ps 1440 MeV, M E( i 660) = (1630 - 1690) MeV ps 1660 MeV, M A(1600) = (1560 - 
1700) MeV « 1600 MeV, and M x(1835) = 1833.7 MeV from the Review of Particle Physics 
[20J. The first radial excited state of the E has not been observed yet, we take the ap- 
proximation M u , = M 3 + Ms(1660) ~ M * +M ^ zjW^Ogoo) -Ma ^ M ^ + 486Me V. The 
strength A and the distribution width /j, are free parameters, we take the values A = 215 
and [i = 200 MeV for the pp bound state as in Ref.|13j. and take the simple replacements 

j^2 j^2 ^2 

MjoS'j Pjfl an< ^ to take into account the flavor SU(3) breaking effects for the 
EE (££(1660)), EE (EE'') and AA (AA(1600)) bound states, respectively. 

We solve the Bethe-Salpeter equations as an eigen-problem numerically by direct it- 
erations, and observe that the convergent behaviors are very good. For the pp, ££, EE, 
piV(1440), IE(1660) and EE' bound states, there exists a solution with A(P 2 = -M x ) = 1 
and Ex < 0. On the other hand, we cannot obtain a solution to satisfy the condition 
X(P 2 = —M x ) = 1 for the AA and AA(1600) bound states. Experimentally, there are AA 
near threshold enhancements in the decays B + -> AAK + , B° -> AAK°,AAK*° [2"T| [22] 
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PP 






piV~(1440) 


££(1660) 




M x [MeV] 


1833.7 


2317.8 


2612.4 


2344.0/2374.0* 


2798.4/2828.4* 


3103.4 


E x [MeY] 


-42.9 


-61.0 


-31.0 


-34.3 


-51.0 


-26.0 


expt[MeV] 


71836.5 


7 2376.3 


7 


7 2376.3 


? 
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Table 1: The masses Mx and bound energies Ex of the baryon-antibaryon bound states, 
the * denotes that the upper bounds of the masses of the iV(1440) and £(1660) baryons 
are taken. 





piV(1440) 


££(1660) 




M x [MeV] = 200MeV,r = 0.55) 


2378.2/2408.2* 






M x [MeV] = 200MeV,r = 0.61) 


2375.3/2405.3* 


2849.3/2879.3* 




M x [MeV] {fx = 200 MeV, r = 0.78) 


2363.2/2393.2* 


2832.9/2862.9* 


3129.3 


M x [MeV] {fx = 200 MeV, r = 0.90) 


2353.2/2383.2* 


2813.9/2843.9* 


3118.1 


M x [MeV] {fx = 200 MeV, r = 1.00) 


2344.0/2374.0* 


2798.4/2828.4* 


3103.4 


M x [MeV] {n = 400 MeV, r = 0.44) 


2378.2/2408.2* 






M x [MeV] (n = 400 MeV, r = 0.47) 


2373.2/2403.2* 


2849.3/2879.3* 




M x [MeV] (p = 400 MeV, r = 0.52) 


2361.8/2391.8* 


2832.3/2862.3* 


3129.3 


M x [MeV] (/x = 400 MeV, r = 0.60) 


2336.8/2366.8* 


2784.4/2814.4* 


3084.4 


expt[MeV] 


7 2376.3 


? 


? 



Table 2: The masses Mx of the baryon-antibaryon bound states with variations of the 
input parameters, the * denotes that the upper bounds of the masses of the iV(1440) and 
£(1660) baryons are taken. 

and above threshold enhancements in the decays B + — > AAK + [23] from the Belle col- 
laboration, the above threshold enhancements can be identified as the J/ip and rj c mesons 
respectively, the decays J/tjj — > AA and ry c — > AA are observed. The AA near threshold 
enhancement may be a AA baryonium state or just a final-state re-scattering effect, more 
experimental data are still needed to identify it. We can study the baryon-antibaryon 
scattering amplitudes in unitary Chiral perturbation theory by taking into account the 
intermediate multichannel baryon-loops (for example, the pp, AA, EE, etc.), and adjust 
the parameters in the phenomenological lagrangian to reproduce the pp baryonium state 
X(1835), and explore whether or not there exists a pole related with the AA baryonium 
state. 

The numerical results for the Bethe-Salpeter amplitudes are shown in Fig.l and the 
values of the bound states are presented in Table 1. From the Table, we can see that the 
new resonances X(1835) and V(2370) can be tentatively identified as the pp and ^^(1440) 
(or A r (1440)p) bound states, respectively, while the mass of the SS bound state disfavors 
identifying it as the X(2120) or JT(2370), because the energy gaps M s g — Mv(2i20) = 
195.4 MeV and M sg - M x(2370) = -58.5 MeV. 

In Ref.[7], Hao, Qiao and Zhang study the h three-gluon glueball with the QCD 
sum rules, and observe that its mass lies in the region of (1.9 — 2.7) GeV, while the 
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Figure 1: The Bethe-Salpeter amplitudes of the baryon-antibaryon bound states where 
p2 = -M%, A, B, C, D, E and G denote the pp, ££, SS, pN(U40), SS(1660) and 
HH', respectively. The line-shapes of the D, E and G are plotted with the replacement 
q 2 //j? — > (q 2 — 10/i 2 )//i 2 to avoid overlapping with the A, B and C. 




Figure 2: The Regge trajectory for the rj mesons. 
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quenched lattice QCD calculations indicate that the pseudoscalar glueballs have masses 
about 2590(40) (130) MeV [21] or 2560(35) (120) MeV [25] . We should bear in mind that it is 
the quenched lattice QCD not the full lattice QCD, as the fermion determinant is neglected. 
In fact, the predictions vary with the theoretical approaches, for example, the mass of the 
low-lying pseudoscalar glueball G is about 2.22 GeV from the quenched QCD Hamiltonian 
in the Coulomb gauge |26j . 2.62 GeV from the constituent gluons model with the Cornell 
potential within the helicity formalism [27], 2.28 GeV from the string Hamiltonian derived 
from the vacuum correlator method [28], 2.19 GeV from the refined Gribov-Zwanziger 
version of the Landau gauge [29], etc. In Ref . [30] . H. Y. Cheng deduces the mass of the 
pseudoscalar glueball G from an rj-rj'-G mixing formalism based on the anomalous Ward 
identity for transition matrix elements, and find a solution Mq = (1.4±0.1) GeV with the 
inputs from the KLOE experiment. If the r/(1405) is confirmed as a pseudoscalar glueball 
one day, it is not necessary that the low-lying pseudoscalar glueballs should have masses 
about 2.5 GeV. The new resonance V(2120) may be a pseudoscalar glueball, and there 
may be mixings among the baryonium states and the glueballs, as the baryon-antibaryon 
bound states can annihilate into three gluons. 

In calculations, we have used a universal coupling constant g 2 (g = g'). In fact, no expe- 
riential data warrants that the coupling constants gNN-w, 9nn(\4AQ)-k an d 5iv(i44o)Af(i440)7r 
have equal values. The experimental value gjv ^^ 40)7r ~ 0.39 — 0.55 [31] and the the- 
oretical value from a special non-relativistic quark model 9iV ^ 1 ^ 40) ' r ~ 0.33 [32] devi- 
ate from 1 obviously, we expect that the values 15^(1440)^(1440)^ I < IffAWTrl- In phe- 
nomenological applications, we often introduce the monopole (or dipole) form-factors 
[33] and the exponential form-factors [M] to parameterize the off-shell effects, and there 
are some form-factors associate with the coupling constants. In this article, we use the 

modified Gaussian distribution g 2 (k) = ^4^^-^ exp ^— j^j , which is assumed to take 
into account the form-factors effectively, and introduce a parameter r (with the value 
< r < 1) to parameterize the difference between the gNN-w and the g./v(i44o)Ar(i44o)7r i 

i.e. g 2 (k) g(k)g'(k) = tA {jp) ex P ( — J?)- We solve the Bethe-Salpeter equations 
with variations of the parameters r and /j, where the flavor SU(3) breaking effects for 
the ££(1660) and EE' bound states are taken into account by the simple replacements 

M 2 A/2 

H fJ-wt anQl respectively, the eigenvalues are presented in Table 2. In calculations, 

we observe that in some regions, there indeed exist solutions in the channels piV(1440), 
££(1660) and HE', and some eigenvalues of the p7V(1440) (or piV(1440)) bound state with 
T = 0.44 - 1.00 and [i = (200 - 400) MeV are consistent with the mass of the V(2370). 

The radiative decays of the J/tp are generally believed to be glue-rich, which can 
explain the branching ratio of the decay J/ip — > 77/ is large (about (5.28 ± 0.15) x 10 -3 ), 
while the branching ratio of the J/ip —?- 77/ is small (about (1.104±0.034) x 10 -3 ) [2D] . The 
observation of the X(1835), V(2120) and X(2370) in the r/' channel not in the 77 channel 
maybe due to the intermediate virtual gluons are flavor-neutral and the r/ meson is mainly 
an SU(3) flavor singlet and has considerable gluon constituent via the axial anomaly. It 
is natural to assume the V(1835) and X(2370) have some gluon constituents, which play 
an important role in the decays to 7r + 7r~r/. 

The hadronic molecular state A which consists of a meson pair or a baryon pair B + C 
can decay through two typical routines, the first one is A— > B + C — > B + E + F + ■ ■ ■ , 
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and the second one is A —> B + C — > E + F + G + ■ ■ ■ , then the decay widths are 
determined by the intermediate process C — > E + F or the annihilation of the B + C, we 
can estimate the widths of the molecular states via the decay mechanisms. For example, 
in Ref.[35], Guo, Hanhart and Meissner take the K(4660) as a i///o(980) molecular state 
considering the nominal threshold of the ip' — /o(980) system is about 4666 ± lOMeV 
[20] . The y(4660) decays dominantly via the decay of the scalar meson /q(980), i.e. 
Y(4660) -)■ V>7o(980) -> ^'tttt, 7//iVKT, and the width of the y(4660) originates from the 
decay of the / (980) mainly. On the other hand, if we take the X(3872) as the D*D±D*D 
molecular state, the decay X(3872) — > D°D°-k (> can occur through the decays D* — > Dir° 
and D* — > Dtt° [36], the narrow widths of the D* and D* mesons warrant that the 
width of the X(3872) is not broad, furthermore, the decay is suppressed kinematically 
in the phase-space. In the present case, the thresholds 2M P = 1876 MeV > 1835 MeV 
and M p + M N{U4Q) = 2408 MeV > 2376 MeV, the decays X(1835) -> pp -> rf V+tT 
and X(2370) -4 pN (1440) , N (1440) p -4 T/vr+Tr^ can take place via the Okubo-Zweig- 
Iizuka super-allowed fall apart mechanism with re-arrangement in the color space. The 
decays V(1835) -4 pp and V(2370) -4 pN (1U0) , N (U40)p occur through the higher 
tails of the mass distributions, and the widths may be large, although the decays are 
suppressed kinematically at the lower tails of the mass distributions. We can search for the 
pN (1440), N(lA40)p enhancements in the radiative decays J/ip — > 7/>iV(1440), 77V(1440)p. 
The recent BESIII data indicates that the Jf (1835) has the width V = (190 ± 9±H) MeV 
[31[4], it is too large for a pure molecular state. A larger mass glueball constituent G besides 
the pp component in the V(1835) is needed to take into account the experimental data, 
the decays G — > ggg — > qqqqqq can take place easily if kinematically allowed, furthermore, 
such glue-rich processes prefer the final-state r]'ir + Tr~ . The conventional qq components 
in the Z(1835) and X(2370) can also lead to the decays X(1835) -4 pp and V(2370) -4 
pN (1440) , N (1440)p with the creation of additional two qq pairs from the QCD vacuum, 
however, the final-state r]ir + ir~ (rather than r]'ir + ir~) is preferred as such processes are 
not glue-rich. It is difficult to calculate the decay widths of the V(1835) and X(2370) 
quantitatively in the framework of the Bethe-Salpeter equation. 

If those bound states presented in Table 1 exist indeed, they can be produced in the 
radiative J '/ip decays, i.e. J/ip — > jgg, gg + qq — ^ pp, EE, SH, piV(1440), EE(1660), those 
bound states can decay to the rjiriT, tjKK, r/'ir-K, r/'KK, rj'rjrj, rfr/rj, ryqq final states. We 
can search for those bound states in the rjirir, r/KK, rj-Kix, r]'KK, rj'rjrf, rfrfr], r/r/rj invariant 
mass distributions in the radiative decays of the J/ip at the BESIII [37J or the charmless 
S-decays at the KEK-B. 

We do not exclude the canonical explanations, the X(1835), X(2120) and X(2370) 
may be the conventional mesons which originate from the confining QCD forces and 
consist of the constituent quark-antiquark pairs with (or without) some gluon compo- 
nents. In Ref.pS], Huang and Zhu take the X(1835) as the second radial excited state 
of the rj', the ground state nonet pseudoscalar mesons are {ir, K, rj, rf}, the first ra- 
dial excited states are {7r(1300), K (1460), 77(1295), r/(1475)}, and the second radial ex- 
cited states are {tt(1800), if (1830), r?(1760), X(1835)}. In Ref.[9], Klempt and Zaitsev 
perform detailed analysis of the properties of the 7/(1295), 7/(1405) and 7/(1475), and 
draw the conclusion that there maybe only one rj state, the 7/(1440), which has mass 
about (1200 — 1500) MeV, and identify the X(1835) as the first radial excited state of 
the rj' . In Fig. 2, we plot the (n,M 2 ) for the 7]' mesons, where the n denotes the ra- 
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dial quantum numbers, the Regge trajectory favors identifying the 7/(1475), X(1835), 
X(2120) and AT(2370) as the radial excited rf mesons with n = 2, 3, 4 and 5, respec- 
tively. The Regge trajectory alone cannot result in definite identification. The decays 
X(1835), X(2120), X(2370) — > r]'Tr + ir~ take place through the emission of a pair of 5-wave 
7r mesons, while the decays X(1835), X(2120), AT(2370) — > r]irir have not been observed 
experimentally yet. Whether or not there exist those decay modes is of great importance, 
further experiments are needed to prove or exclude the possibility. 

4 Conclusion 

In this article, we study the mass spectrum of the baryon-antibaryon bound states pp, £E, 
EH, AA, piV(1440), ££(1660), EE' and AA(1600) in the framework of the Bethe-Salpeter 
equation with a phenomenological potential. The numerical results indicate that the pp, 
££, EE, pA^(1440), ££(1660), EE' bound states maybe exist, and the new resonances 
X(1835) and X(2370) can be tentatively identified as the pp and pA^(1440) bound states 
respectively with some gluon constituents, while the new resonance AT(2120) may be a 
pseudoscalar glueball. The other bound states predicted in this work may be observed 
experimentally in the future in the radiative decays of the J/ip at the BESIII or the charm- 
less -B-decays at the KEK-B. On the other hand, the Regge trajectory favors identifying 
the r/(1475), X(1835), X(2120) and X(2370) as the excited rf mesons with n = 2, 3, 4 
and 5, respectively. 
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